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Ce cours de mécanique des milieux continus est à la base de l'enseignement de mécanique à l'ISITV. Les 
notions abordées ici, transport de champs, lois de conservation, ..., seront reprises ultérieurement en 
mécanique des solides et mécanique des fluides. Dans une première partie, nous aborderons les notations 
tensorielles et vectorielles indispensables à toute étude scientifique, puis dans une deuxième partie, nous 
étudierons la cinématique des milieux continus. Après avoir introduit la modélisation des efforts et les lois de 
conservation par le principe des puissances virtuelles, nous appliquerons ces lois de conservation aux lois de 
comportement de l'élasticité linéaire (en mécanique des solides) et aux lois de comportement des fluides 
newtoniens (en mécanique des fluides). 
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Notations tensorielles 


NOTATIONS TENSORIELLES 


1 Vecteurs et tenseurs 

Avertissement: L'objectif de ce chapitre, est de familiariser les étudiants avec les notations tensorielles. Afin 
d'en simplifier le contenu, nous ne considérerons que des bases orthonormées. 

1.1 Notations 


1.1.1 Vecteur 

Dans un espace euclidien £ à trois dimensions, soit e v e 2 , e g une base orthonormée. Un vecteur V est 


représenté par ses composantes V v V 2 , V 3 

v = vA + vA + vA = èvë t 


2=1 


( 1 . 1 ) 


En utilisant la convention de sommation, ou convention d'Einstein, on écrit 

V = Ve. 

l l 


( 1 . 2 ) 


où, chaque fois qu'un indice est répété, il convient de faire varier cet indice de 1 à 3 et de faire la somme. Dans 
l'expression (2) l'indice i est un "indice muet". 


En notation matricielle on écrira parfois 


V = 



A 

K 


(1.3) 


et le vecteur transposé 

f r ={f}'=(f)=(v;Kn) 


(1.4) 


1.1.2 Application linéaire de £ dans £ 

Soit A une application linéaire, dans la base e ] , e , e 3 . Cette application est représentée par une matrice 3x3 
notée fyl]: 

Ai A* As 
Ai A 2 A 3 
Ai A 2 As. 

Si W est un vecteur tel que W = AV , alors les composantes de W sont données par 

w i = A Ai + AsA + As A 
*A = AAi+AA;+AsA 
^3 = AAi + A^ + AsA 

et en utilisant les conventions de sommation où j est un indice muet 
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W = AV 

• y j 

et en notation vectorielle 

M = MW 

On définit les symboles de Kronecker par 

1 si i = j 
0 si i ^ j 



( 1 . 5 ) 


( 1 . 6 ) 


En particulier l'application identité 1 est représentée par la matrice 


«u 


6 * 


1 

0 

O . 




= 

0 

1 

0 





0 

0 

1 


La composition de deux applications linéaires se traduit par le produit de leur matrice représentative, c'est-à- 
dire 


C — A o B ou encore [(?] = [.4] [i?] 


et en notation indicielle 


C = A.B 

IJ IK KJ 


( 1 . 7 ) 


1.1.3 Formes bilinéaires 

Soit A une forme bilinéaire sur £, c'est-à-dire une application bilinéaire de £xÇ dans K. Dans la base 
e v e,, e 3 elle est représentée par une matrice A telle que 


a{v,W^ 


A.V.W 

V * J 


( 1 . 8 ) 


ou en notation matricielle 

A^VW^fr^w} 

En particulier, la forme bilinéaire représentée dans toute base par les symboles de Kronecker est le produit 
scalaire. Si ( e v e 2 , e 3 ) est une base orthonormée, alors 


e. • 

l 


e. 

j 



et le produit scalaire de deux vecteurs est donné par 

V • W = Ve. • W.e. = V.W e.-e.= 6..V.W. = VW 

ii jj i j i j ij i j i i 


ou en notation matricielle 

v -w = ly){w) 


1.1.4 Tenseurs 

1.1. 4.1 Tenseur du second ordre 

Un tenseur du second ordre T est un opérateur linéaire qui fait correspondre à tout vecteur V de l'espace 
euclidien un vecteur W de ce même espace. 
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W = T[V 


Cet opérateur peut être représenté par une matrice 3x3, notée T ou 


ou T , telle que 


W. = TV 

* tj J 


ou en notation matricielle 


W 


V 


ou 


W = TV 


* Un tenseur est dit symétrique si T. = T. 

IJ Jl 

* Un tenseur est dit antisymétrique si T = —T 

' n y jt 

* Un tenseur est dit isotrope si T — t6.. 

y y 


* On peut toujours décomposer un tenseur en une partie symétrique et antisymétrique 


T = T +T 


rj - r rj-\S rpA 

ou ÿ ÿ ÿ 


T s =-(t.+t) T a =-[T -T 


v n \ v ji f tj 

avec z et 


1.1. 4. 2 Tenseur d'ordre supérieur 

On peut définir un vecteur V par ses composantes V , ou par les coefficients de la forme linéaire 

X—>X-V = Xy., car la base choisie est orthonormée (voir les notions de vecteurs covariants et 
contravariants). 

On peut alors considérer le vecteur comme un tenseur du premier ordre. 

De même, une fonction scalaire peut être considérée comme un tenseur d'ordre zéro. 

Un tenseur du troisième ordre S est un opérateur linéaire qui, à tout vecteur Z fait correspondre un tenseur 
du second ordre T . 

T = S ( Z ) ou encore T.. — S. jk Z k 


1.1. 4. 3 Produit tensoriel 

On définit le produit tensoriel du vecteur U par le vecteur V , noté U ®V , comme le tenseur d'ordre deux, 
défini par la forme bilinéaire qui aux vecteurs X et Y fait correspondre [u-x\[v-y\ 

Les 9 produits tensoriels e. ® e définissent une base de l'espace vectoriel des tenseurs d'ordre deux, si bien 
que l'on peut écrire un tenseur T comme 

T = Te ®e 

y * j 

ou encore, par exemple, 
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u ® v = u.v.e. <8> e. 

‘ i i j 

1.1. 4. 4 Contraction et produit contracté 

Soit le produit tensoriel A® B ® C , on appelle contraction, l'opération qui lui fait correspondre le vecteur 

A(B- C) . Le produit contracté d'un tenseur d'ordre 4 R et d'un tenseur d'ordre 3 S est défini par le tenseur 
d'ordre 5 


U , V, 

uv n 

uv n 

1 1 

1 2 

1 3 

u n v. 

u n v n 

UV^ 

2 1 

2 2 

2 3 

u n v. 

u n v n 

uv„ 

3 1 

3 2 

3 3 


R- S = \R m e. 


S e 

pqr p 


= R S e. 

ijkm mqr i 


Le produit doublement contracté d'un tenseur d'ordre 4 fi et d'un tenseur d'ordre 3 S est défini par le 
tenseur d'ordre 3 

R : S = (üL,e. ® e. ® e. ® e.) : (s e ® e <8> e ) = R S e. ® e. ® e 

y ijkl i j k l ) \ pqr p q r ) ijnm mnr i j r 

Par exemple, le produit doublement contracté de deux tenseurs d'ordre 2 T et T' est le scalaire 

T \ T' — (Te. <g>e.) : (t'J ®e ) = T T'. 

\ v * j j y- 1 pi p « ) vp 

1.2 Changement de repère 


1.2.1 Matrice de passage 

Soit e , e 2 , e 3 une base orthonormée et e 1 , e', e' une autre base orthonormée. 


On définit la matrice de passage Q telle que: 


A ^11^1 + *3l2®2 + ^13^3 
Ê 2 _ ^2i e i Q'22 e 2 Q>3 Ê 3 

K = QA + QA + QJ 3 

ou encore, en notations indicielles 

e' = Q e 

l 'v y j 

et en notation matricielle 

{e'} = [Q] {e} 

Les deux bases étant orthonormées, on doit avoir 

Ai = A ■ A' = QiïAk ' QjÂ = = * 3 * * 3 ,* 

ce qui montre que la matrice inverse de Q est Q 1 . En particulier on tire la relation inverse: 

e. = Q e' 


1.2.2 Vecteurs 

Soit V un vecteur de composantes V dans la base e,, e„, e, et V' dans la base e', e', e( . 

1 i L ' 2 ' o l 1 ' 2 ' à 

v = v.ë. = v.ë' 

ii ii 
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En utilisant la matrice de passage 

V = VÏ=V ! Q h l k 


soit 


v; = va et V=vü 


ou encore, en notation matricielle 



Remarque: le produit scalaire est un invariant, c'est à dire que cette fonction est indépendante du repère 
choisi. 

En notation indicielle 

V'w' = v\v! = VQ.W.Q. = S VW. = VW = VW 

k k i ^ ki j ^ kj i] i j il 

et en notation matricielle 


V'W = (V')\W'\ = 


Q\\W 


_[«]{? 

= (^)M t [q]{^} = (^){^} = v- w 


1.2.3 Application linéaire 

Soit A une application linéaire, de composantes A dans la base ëj, e 2 , e 3 . et A', dans la base e', e', e' . 
En notation indicielle 

rç'= A X = Q = QAm v m = QAmQ^y: 

d'où 

A' =Q.A. Q. 

ik ^ ij jm ^ km 

et en notation matricielle 

F'} = Mp) = MH = MMp} = MMM'H 

soit 

M=MM|o] r 


1.2.4 Forme bilinéaire 

Soit A une application linéaire, de composantes A dans la base e v e 2 , e 3 . et A!, dans la base e', e', e’ 3 . 


A(V, W) = A.VW = A'.VW = A Q.V'Q W' 

\ J / ij i j ij i j ij ^ la k ^mj m 

soit 

A! = A Q Q . 

km ij ^ ki ^mj 

et en notation matricielle 

A{V,W) = /v\[A]p} = = 

[«r{pfMi«rp}=F)NHi«rF} 
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1.2.5 Tenseur d'ordre 2 

Soit T un tenseur d'ordre 2, en notation indicielle 


T = Te 

IJ l 


= Tl'. < 

y * 


e'. = T Q e' <8>Q e' = T Q Q ~e[ ®e 

i y w Kl K 'mi m ».i w lei Wi mn le i 


ij^ki^mj k 


puis 


T' = TQQ 

km m ^ ^ m 


2 Permutations et déterminants 

2.1 Les symboles de permutation 

On introduit les symboles de permutation 


£ 


ijk 


+ 1 si i,j,k est une permutation paire de 1,2,3 
— 1 si i,j,k est une permutation impaire de 1,2,3 
0 si deux indices sont répétés 


Ces symboles représentent le produit mixte des vecteurs de base 


"ijk 


e.. e.. e, 

J ’ k 


e... sont les composantes d'un tenseur du troisième ordre, qui représente, par exemple, la forme trilinéaire 

IJK 


produit mixte: 


u,v,w\ = e Vs uyw i 


Avec un peu de patience on peut démontrer les résultats suivants 


= Det 

ijk irnn 


vm m 


6 .. 8 

il i 

8 , 8 . 8 . 

jl jm jn 

^kl ^km ^kn 


£...£. = 8 . 8 . — 8 . 8 

iik imn im kn m km 

£ ijk £ ijk~~ ^ 


2.2 Déterminant d'une matrice 


Les symboles de permutation permettent le calcul du déterminant d'une matrice par 

£...Det(A) = £ A. AA 

ijk v ' mnp im jn kp 


(1.9) 


ou encore 


Det(A) — — £...£ A A A, 

' ' r* ijk mnp im jn kp 


On peut également déterminer l'inverse d'une matrice 
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B = A~ 2 et B = 

]t 2Det(A) 


2.3 Polynôme caractéristique 


£ . £ . A A 

tmn jpq mp nq 


Les valeurs propres d'un tenseur du second ordre sont obtenues par la résolution de l'équation caractéristique 

P( A) = Det(A - XI) 

soit en développant 

— £..,£ (A. - XS VA. -XS. )(A — XS.) = 0 

n nk mnv ' vm im ' \ m m ' x ko ko ' 


6 


ou encore 


P( X) = I 3 - X I 2 + A 2 /j - A 3 


avec 


I = —£...£ A. A A = Det(A) 

d g tjk mnp vm jn kp v ' 

I = -( A A. - A..A.) = -\(Tr Af -Tr A 2 1 

2 c\ { il JJ ij Ji) n [' 7 J 

I = A — Tr A 

1 u 


I v I 2 , 1 3 sont appelés les invariants fondamentaux du tenseur A. 

2.4 Adjoint d'un tenseur antisymétrique 


Soit 11 un tenseur antisymétrique 


fl = 


0 

-fl. 


^12 ^31 


fl 


fl -11 

31 23 


on peut également lui associer le vecteur 


U! = 


ÜJ, 


fi,. 

1 


23 

ÜJ n 

— 

Q 

2 


31 

LU, 


f2 

3 


12 


soit 


fl = 


Le vecteur lu est le vecteur adjoint du tenseur antisymétrique fl . En notation indicielle on a: 


fl = E...U), 

ij ijk k 

ÜJ =->-£. fl 

i 2 y* jk 


( 1 . 10 ) 
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3 Calcul vectoriel et analyse vectorielle 

3.1 Calcul vectoriel 

Le produit vectoriel 

c = a A b 


s'écrit en notation indicielle 


ce. = £.,a.b.e. 

il ijk j k, i 


On peut montrer que 

\a Ab\Ac = \a- c\b — \ b- c\ a 


a A 6 • c A rf = [a- c)\b- d\ — [a- d\[b- c 


3.2 Analyse vectorielle 

On note d'une virgule la dérivée partielle, soit 

exprimés dans un repère cartésien orthonormé. 



Les opérateurs exposés dans cette partie seront 


* Soit / une fonction scalaire 

Le gradient d'une fonction scalaire est un vecteur 


grad f = Vf = f. e. 


df_ 

dx 1 

df. 

dx 2 

df 


dx , 


Le laplacien d'une fonction scalaire est un scalaire 


d 2 f d 2 f 
A f = f..=fLL + fLL. 

■" dx: dx l 


5V 

dx\ 


* Soit v un vecteur 

La divergence d'un vecteur est un scalaire 


Div v = v. 

l,l 


dv l 

dx 1 



dx 3 


Le rotationnel d'un vecteur est un vecteur 


rot v = V A v = £ ... v, e. = 

IJK K, J l 


dx 2 

dv 1 

dx 3 

dv 2 

dx 1 


dV 2 
dx 3 

OUa 

dx 1 

dv 1 

dx 2 


Le gradient d'un vecteur est une matrice 
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V v = v e. 


dv 1 

dv 1 

dv 1 

dx 1 

dx 2 

dx 3 

dv 2 

dv 2 

dv 2 

dx 1 

dx 2 

dx 3 


dv s 

dV 3 

dx 1 

dx 2 

dX 3. 


Le laplacien d'un vecteur est un vecteur 


A v = v. . e. 

i,]j t 


' d\ 


d 2 Vl 


d\ 

dx\ 


dxl 



d\ 


d\ 


d\ 

dx\ 


dxl 


dx 2 3 

d\ 


d\ 


d\ 

dx 2 

~r 

dx\ 

~r 

dx l 


* Soit T un tenseur du second ordre 

La divergence d'un tenseur est un vecteur 


Div T = T . . e . 

>ij « 


* Quelques formules utiles 

Div^f aj = f Diva + a ■ grad f 
Div \a/\b) = b-rota — a-rotb 


Div rota = 0 


rot grad / = 0 


gradlfg) = f grad g + g grad f 


rot \fa\ = f rot a + grad f A a 
Div {grad /j = A / 
rot^rotaj = gradyDiv aj — A a 

3.3 Transformation d'intégrales 


Av o _ 

Au„ 


dT n 


dT u 


vt 13 

dx 1 


dx 2 


dx 3 

dT 21 


dT 22 


dT 

23 

dx 1 


dx 2 


dx 3 

VT, 


vt, 2 


dT 

33 

dx 1 


dX 2 


dX 3 . 


Soit fl un domaine borné et dfl sa frontière, de normale n 
Soit (f> une fonction scalaire, alors 


JJaat ndS = IIfn9 rad <t> dV 

Soit A un vecteur, alors 

ff gn A.n d S = fff Si mv(A)dV 
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Soit T un tenseur, alors 

JJ dQ T.ndS = JJJ n Div(T)dV 

Soit dfl un domaine plan de normale n , de frontière T . Soit U un vecteur défini sur ce domaine. Si r est le 
vecteur unitaire tangent à T , alors 

fLn rot (U) ■ 11 dS = f r U ■ T dl 

Tous ces résultats sont issus du théorème de la divergence 

ffaa^i dS = UfJm dV 


4 Formules essentielles en Mécanique des Milieux Continus 

4.1 Coordonnées cartésiennes orthonormées 

OM = xe +ye + ze 

x a y z 


: Soit v = v e + v e +ve un vecteur, alors 

xx y y z z 


— d y _ _ 

Vm = Vu = — - e. (g) e. = v. . e. ® e. = 

v ' Q x • 1 hj • 3 


dv 

x_ 

dv 

x_ 

dv 

X 

dx 

dy 

dz 

dv 

y_ 

dv 

y_ 

dv 

y_ 

dx 

dy 

dz 

dv 

dv~ 

dv . 

dx 

dy 

dz 


et 


dv. 


dx. 


divv = — = v. . = Tr\ grad{v) = Vu : / = — — H — + 


dv dv_ dv 


dx dy dz 


Av = div V(u) = 


d 2 v. 

l 

dx.dx. 

J J 


e = v. .. e = Av e + Av e + Av e. 


* Soit f une fonction scalaire, alors 


df 


grad(f) = V/ = — e. = f. e. = 


et 


* n . d 2 f , 


/ dx.dx . 

3 3 


* Soit T = T., e. <g> e. = 

y * 3 


T T T 

xx xy xz 

T T T un tenseur symétrique du deuxième ordre, alors: 

yx yy yz\ 

T T T \ 

zx zy zz 
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= dT 

div(T) = — ^ e=T 
dx. ' 


e. = 


9T 

XX 


dT 

xy 


dT 

xz 

dx 


dy 


dz 

dT 

yx 


dT 

yy 


dT 

yz 

dx 


dy 


dz 

dT 

zx 


dT 

zy 


dT 

ZZ 

dx 


dy 


dz 


et 


= d 2 T. 

AT = lJ - 


dx k dx k 


e =T. e. 

J ij,kk i 


> e. = 
j 


AT AT AT 

xx xy xz 

AT AT AT 

yx yy yz 

AT AT AT 


4.2 Coordonnées cylindriques 


7-7 dOM _ 

(JM = re + ze et = e 

r z t-\ r 

or 


d(OM ) = e dr + rd0e fl + ? d^ 


1 90M 
r d 6 


dOM 

dz 


de 

r 

— n 

03 

7 * 

0 

—P 

de 

z 

— D 

dr 

de 

r 

— U 
— P 

1 ^ 

Ü*iQ3 

"S 

' 1 

dr 

de 

Z 

— U 

— D 

de 

e o ' 

’ de 

C 

r 

’ de 


de 

r 

= 0 . 

I ^ 

II 

0 

de 

= 0 

dz 


dz 


’ dz 



: Soit 1; = v e + v„e„ +ve un vecteur, alors 

r r 0 0 z z 


grad(v ) = Vu = 


1 

9r r 


dv„ 


1 


9r r 

dv 

dr 


dv 

~ V r> 

de e 

dv e 

h V 

de 

1 dv, 

r de 


dv 

r_ 

dz 

dz 
dv „ 
9z 


et 

dit; u = Tr |v(u)J 

Au = dùufvuj = 

* Soit / une fonction scalaire, alors 


v dv l 9u„ 9u_ 


1 = 

r 

+ -^ + - 
9r r 

0 

de 

_| — 

dz 

2 dv e 

V 

r 

~p -4— 

Av g 

2 dv 
_l_ A r_ 

r 2 de 

r 2 

r 

r 2 90 


e, + Av e 


, 9/_ 1 9/ _ 9/_ 

grad(f) — Vf — — e 4 e. + — e 

ar r a# az 


et 


A/ = diu(v/) 


9 2 / ^9/ ^97 97 

9r 2 r 9r r 2 90 2 97 
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! Soit T = 


T T „ T 

it rQ rz 

T T T 

Or 60 Oz 

T T „ T 


un tenseur symétrique du deuxième ordre, alors: 


div(T) = 


dT 1 d 7 


d i 


de 


dT 

rz_ 

dz 


T — T.. 


dT 1 d T dT 


dr 

dT 


r de 
1 dT 


de 


dz 

dT 

22 

dz 


r 

r 

T 


4.3 Coordonnées sphériques 


7777 _ , dOM _ 

OM = re et = e 


dr 

de 


= 0 


de 


de 


1 dOM 


1 dOM 


dr r r de 
d(OM ) = edr + rdde g + rsind d<\> e ' 
de 


rsind d(j> 


d<j> 


■ = sinde 


de, 

dr 

de g _ g 

de ~ e " 

de, 

— - = cosde, 
d4> * 


de, 
<P_ 

dr 

de, 


= 0 


de, 
0 

d(j) 


= 0 

de 

= sinde — cosde. 


Soit v = v e + v,e, + v,e, un vecteur, alors 

r r 6 9 <p <p 


grad[v) = V« = 


dv 1 

r_ _ 

dr r 

I 

dr r 

d y à 

dr 


et 


- - ~ dv 

divv = Vv : I = — - 
dr 


dv 

r_ 

de 

dv 

— + 

de 

1 dv t 

r dd 


1 dv 


sine d(j) 
1 dv, 


sine d(j> 

1 dv é 


siriO d(j) 


— cotgdv 1 
cotgdv + v 


v 1 dv, 
■2 — H fl 


r r dd rsind d(j> 


i dv v 

* + cot gd — 
r 


Av = div\X7(v) 


Av — — 


Av e + — 
r 


1 < 9 (sin dv) 1 dv 


sin d dd 


sin 0 d(f> 


dv 


cos 


d dv. 


dd 2 sin 2 d sin 2 d d<p 


Av 


2 dv n dv. v 

• -i + cotgd — - — 

* r 2 sin d d(f> dcj) 2 sin d 
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: Soit / une fonction scalaire, alors 


grad(f) = 


df_ 

dr 

1 df_ 
r d9 
1 df 


rsin# d(j) 


et 


A , d 2 f 2 97 1 .df 1 d-f 

Af = div ( grad(f) \=—±. + --±. + — cotgd + , — 

I dr" r DO r 30 r" sin 9 d(f) 


: Soit T 


/ / / 

rr rO rcf) 

T T T 

Or 00 0 <p 

T T T 

(j>r cj)0 (jxj) 


un tenseur symétrique du deuxième ordre, alors: 


div(T) = 


dT i dT a 

rr _| 1_ ru 


1 


dT , 


dr r 86 r sin 0 d(f> r 


+ 1 (^- T oo-T H , + T re ™t 9 e) 


dT er , i dT„ 


1 


dT 

ûj 


L + - ( (T m - T. . )cotgd + 3 T, ) 
dr r dd r sin 9 dcj) r ' ^ ' 

dT \ 


dT *r , 1 dT* 


1 


dr r d9 rsm9 dcf> 


+ i (2 + 


4.4 Comment retrouver les formules 


Nous nous plaçons par exemple en coordonnées cylindriques. On note 

v — v e + u.e. +ve= v.e. avec i == r, 2 et , i = — , — 

or r oü oz 

Donc, avec cette convention 


e r,o =— et e g 
r 


e 

r 

r 


Chercher le gradient d'un tenseur consiste à augmenter l'ordre de ce tenseur, soit 


V(**) = (**) j 0 e. 

Si on applique cette remarque à un vecteur, on obtient: 

V(w) = (v.e.). 0 e. 

En n'oubliant pas de dériver les vecteurs de base, car nous sommes dans un système de coordonnées 
cylindrique, 

\7v = v. . e. 0 e. + v. e. . 0 e. = v. . e. 0 e. + e.„ 0 e„ 

= v i,j e 4 ® e . + v r e r _ fl 0 e, + v, 0 e 9 

v - - v _ 

= e ® e + — e ® e - e 0 e 

J rjr» rjr» 

Pour obtenir l'opérateur divergence, il suffit de contracter doublement avec le tenseur unité d'ordre 2, 

div(**) = V(**) : 1 

soit dans le cas d'un vecteur: 
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- - v v dv 1 dv a dv 

div(v) = V« :l = i>. — H H H *- 

w w r r dr r dO dz 


et donc l'opérateur Laplacien pour un scalaire 

A <p = div(yipj = ip ü + 


d 2 ip 1 dtp 1 d 2 <p d 2 <p 
r dr 2 r dr r 2 d9 2 dz 2 


Appliquons maintenant cette méthodologie à un tenseur d'ordre 2. 

V(r) = (rioe.) 


= T e. g) e. g) e. + T., e . , g) e. (g) e, + T. e. g) e. g) e. 

î?,*; i j k ij i,k j k ij i j,k k 

= T., e. g> e. g) e. + T., e. . g> e . g> e. + T. e. g> e. . g) e. 

ij,k i j k ij i,9 j 9 ij i j, 9 9 

T T 

= A ® ~ è ] ® +ye 8 g> e. g>e fl --^-e r g> e. g> e fl 


T. 


W 


— e. g) e. g) e. e. g) e g) e„ 

r * 8 0 r > r 8 


Pour obtenir la trace de ce tenseur d'ordre 3 on contracte les deux derniers indices: 


div 


— T T T. 

V(T) : 1 = T.. . e. H — —e — e H — — e. 

' 2 11,1 % 9 r i 

J J J y» ry* 

dT . 1 321 . OT T 


■H 


+ 


+ 


dr r d6 dz r r 
dT er + 1 ?T„ , 32;, , T g t T gr 


dr r d9 dz 

dT r l 3T„ 321 


_ _j zz_ _|_ zr 

dr r d9 dz r 


On peut donc maintenant retrouver l'opérateur Laplacien d'un vecteur : 

Av = div Vw 


v.„-- 


V a a + ‘ 


- . T.e -> 0,0 - 

= v. ..e. H — e, e 

jrj % 9 r 


-e. -■ 


Av 

~p 1 

> 2 t’„ 

Au 4- ü 

r r 2 d9 r 2 

r 

—a c/ n 

9 r 2 30 r 2 


e H —e. 

r i 

r r 


e a + Av e 

9 z z 
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5 A retenir 


Convention de sommation : 

V = Ve. 

l l 

Produits tensoriels : 

u ® v = u.v£. ®e.— 

< i < i 

Symboles de permutation : 

+1 si i , j, k est une permutation paire de 1 , 2 , 3 
— 1 si i,j,k est une permutation impaire de 1 , 2,3 
0 si deux indices sont répétés 

Produit vectoriel : 

c = a A b = s... a b. e 

tjk j k i 

Quelques opérateurs : 

Div v — v. . , rot v = V A v = £.,v, . e. , V v = v. . e. ® e. , Div T — T e. 

i,i ijk k,j i ' i,j i j ' ij,j i 

En systèmes de coordonnées cylindrique ou sphérique, mieux vaut utiliser un formulaire ! 


£ ijk 


U , V. 

u,v n 

UV n 

1 1 

1 2 

1 3 

UJ), 


UV 

2 1 

2 2 

2 3 

uv. 

uv 

UV 

3 1 

3 2 

3 3 
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1 Le mouvement et ses représentations 

1.1 Configuration 

L'espace physique est rapporté à un repère orthonormé direct (0,e v e 2 ,e ) . L'ensemble des particules ou 
points matériels constituant le milieu continu étudié, occupe à chaque instant t, un ensemble de positions dans 
l'espace: c'est la configuration du système à l'instant t, noté O(f) (d'intérieur fi(£) et de frontière dfl(t) ). 

On introduit aussi la notion de configuration de référence: c'est la configuration particulière du système à un 
instant t fixé. Souvent on prendra O n = 0(0) , et on parlera alors de configuration initiale. 

Toute particule M () de 0 () est repérée par son vecteur position X(l') dans la configuration de référence. Toute 
particule M de 0(t) est repérée par son vecteur position x(t) dans la configuration actuelle (à l'instant t). 


La position de chaque particule M sera donc déterminée si on connaît sa position dans la configuration de 
référence et une fonction $ telle que: 


$ définit le mouvement par rapport à (0,e v e 2 ,e ) ■ On devra donc déterminer trois fonctions scalaires, telles 
que: 


Dire que le milieu est continu, c'est dire que $ est une fonction continue et biunivoque de X . On supposera 
que <[> est différentiable. Le déplacement par rapport à la configuration O , à l'instant t, de la particule M (] est 
le vecteur 



m 


Figure 1 : Configurations de référence et actuelle 



( 2 . 1 ) 



( 2 . 2 ) 


ü(X,t) = x{X,t)-X 


(2.3) 
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1.2 Variables de Lagrange et variables d'Euler 

Une grandeur attachée à une particule (masse volumique, vitesse,...) peut être définie, 

- Soit en fonction de X et t : variables de Lagrange 

- Soit en fonction de x et t : variables d’Euler 

Le vecteur vitesse d'une particule M est défini par 


V(X,t) 


dOM 

dt 


d$(X,t) 

dt 


Le vecteur accélération d'une particule M est défini par 


r(i,f) 


dV(X, t) 
dt 


d 2 $(X,t) 

dt 2 


(2.4) 


(2.5) 


1.2.1 Trajectoire 

On appelle trajectoire d'une particule, la courbe géométrique lieu des positions occupées par cette particule au 
cours du temps. x(t) — est une représentation paramétrée en temps de la trajectoire. Par définition 

de la vitesse, 


V(x, t) 


dOM 

dt 


dx, _ dx„ _ dx , ^ 

— e. H -e H -e 

dt dt dt 


les trajectoires peuvent être obtenues par la résolution des trois équations 


dx. 


dx n 


dx , 


V 1 {x 1 ,x 2 ,x 3 ,t) V 2 (x v x 2 ,x 3 ,t) V 3 (x v x v x 3 ,t) 


— dt 


( 2 . 6 ) 


1.2.2 Lignes de courant 

A un instant donné, on appelle lignes de courant du mouvement, les lignes qui sont en tout point tangentes au 
vecteur vitesse de la particule située en ce point. Soit pour t fixé, deux équations: 

dx 1 dx 2 dx 3 

V 1 (x v x 0 _,x 3 ,t) V 2 (x v x 2 ,x 3 ,t) V 3 (x v x 2 ,x 3 ,t) ( 2 - 7 ) 

Remarque : Pour un mouvement stationnaire (ou permanent) V(x, t) = V(x) . Les lignes de courant et les 
trajectoires sont confondues. 

1.3 Dérivées particulaires 


1.3.1 Définition 

Lorsque l'on suit une particule dans son mouvement, la grandeur A attachée à la particule ne dépend que de 
t. Par définition, on appelle dérivée particulaire de A à l'instant t, la dérivée de A par rapport à la seule 
variable t. 

En variables de Lagrange: A = A(X,t) 


dA 

dt 




dt 


En variables d'Euler: A = A(x, t) 


( 2 . 8 ) 
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, . dA . dA . 

oL4(a;, t J = (i, t)at H (a:, t)dx 

i 

cL4 . cM . cM . dx 

— (z, £) = (æ, f) H (x, t ) — - 

dt dt dx dt 

3 

dA dA . cM . 

— (x, f ) = (x, t) H (a-, t)V 

dt dt 1 dx. K ' 3 


ou encore 


" = 21 + |i /.vu 

dt dt 


1.3.2 Application à l'accélération 

ï ( i, t) = ^m = §HV-V\v 


que l'on peut également écrire 


dV 1 - -* — 

T (x,t) = + -W +rotV A V 

dt 2 


(2.9) 


( 2 . 10 ) 


2 Déformation d'un milieux continu 

2.1 Notion de déformation 

On dira qu'un milieu continu en mouvement subit des déformations si les distances relatives des points 
matériels varient au cours du temps. 

En différenciant (2.1), on obtient: 


dx{t) = dX 


d<S> 

dx.e. = -dX e. 

' ' dX 3 ' 


On note F l'application linéaire qui fait passer de l'espace vectoriel dans lequel peut varier dX dans l'espace 
vectoriel où varie a priori dx. Cette application linéaire, appelée tenseur gradient ou application linéaire 
tangente, permet donc le passage de la configuration Q ( . à la configuration Çl(t) . 



En notation indicielle, 





dx l 

dx 1 

dx f 



dX 1 

dX 2 

dX s 

<9$ 

dx. — 

= — — soit F = 

dx 

dx 0 

dx 

i 




dX. 

dX. 

dX 1 

dX 2 

dX 3 



dx 

dx * 

dx 2 



dX 1 

dx 2 

dx 3 


( 2 . 11 ) 
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2.2 Tenseur des déformations 


2.2.1 Définition 

Le tenseur gradient décrit la transformation locale au voisinage d'une particule donnée. Afin de rendre compte 
des déformations, c'est à dire des changements de forme autour de cette particule, on s'intéresse à l'évolution 
du produit scalaire de deux vecteurs matériels pris respectivement dans les deux configurations f2 n et Çl(t ) . 

Considérons trois particules voisines X , X + dX , X + dX' . Après déformations, elles occupent dans f l(t) 
les positions respectives x , x + dx , x + dx' . 



d'où sa variation autour de la transformation 

dx ■ dx' — dX ■ dX' = 


dx. dx'. 

k k ^ 


dX dX' 


dX.dX' = \f,.F,. -ô. ldX.dX' 


ki kj 


soit 


dx ■ dx' - dX ■ dX' = 2 dX e dX' 


en posant 

= _ 1 
2 

L'application linéaire g est appelée tenseur des déformations. Cette application est symétrique mais dépend 
bien sûr de la base (0,e v e 2 ,e 3 ) initialement choisie. 


F (X,t)F(X,t)-l 


2.2.2 Remarques 

* S'il n'y a pas de déformations, alors e = 0 (et inversement). 

* C = F F est appelé le tenseur des dilatations. Ce tenseur est symétrique. 
On peut démontrer: 

Théorème 1: Les valeurs propres de C sont strictement positives. 


Théorème 2: Det 



>0 


Vf 


Théorème 3: e est symétrique et possède les mêmes vecteurs propres que C . 
* Variation de longueur 

Soit dX' — dX = dl n e et cm = dl, alors 

0 x ' 
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dx ■ dx' - dX ■ dX' = dl 2 - dl 2 = 2 dX e dX' = 2dl:e 


ou encore, si les déformations sont petites 

dl r, — I — „ dl - dl n 


dl 


-»/l -b 2s 1 4" e — 6 

y XX XX Xi 


dl 


0 0 

e ti représente au premier ordre la variation de longueur dans la direction x . 

* Variation d'angle 

Soit dX = dl n e , dX = dl e , alors 

0 x ' 0 y 


dx-dx-dX- dX = cos Bdldl' =2 dX e dX' = 2 dis 

U xy 


ou encore, 


2e = cos 


^ + 2e ™ 


1 + 2e 


donc e^ représente au premier ordre la variation d'angle entre les directions x et y . 

2.2.3 Autre écriture 
D'après (2.3) et (2.1) 


F(X, t ) = — (X, t) = 1 + — (X, t) 

\ t gx K t gx \ 


soit 


- 1 

e = — 
2 


+ d±(X,t)+d±(X,Ax,t) 

dx dx dx dx 


ou encore en notation indicielle 


e.. = — 

11 2 


du du du, du, 

% I J | K K 


dX dX dX dX 

J • > J 


(2.13) 


2.2.4 Cas des petites perturbations 

Cette hypothèse correspond au cas où u(X,t) et 


—(■Ci) 

dX 


sont petits. 


En reprenant (20) et en ne retenant que les termes d'ordre 1, on obtient: 


Eh pp 


%kx,t) + 9u_ {X ,t) 
dx dx 


ou encore en notation indicielle 


£ ijHPP 2 


du. du. 


dX dX 


(2.14) 
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2.3 Conditions de compatibilité 

A tout déplacement Ü on fait correspondre une déformation e . On peut aussi se poser le problème inverse. 
Ce problème est dit 'problème de compatibilité géométrique d'un champ de déformation', ou encore 'problème 
d'intégrabilité d'un champ de déformation'. 

Les conditions de compatibilité peuvent être établies dans le cas général, cependant nous ne les établirons que 
dans le cas des petites perturbations. 


Décomposons maintenant le gradient des déplacements en une partie symétrique e et une partie 
antisymétrique tu . 


- 1 

LJ = — 


- du (X,t) 

1 

LJ.. = — 

du. du. 

2 

dX 

dX 

v 2 

dX dX 

3 1 ) 


On a 


^ij.k ^ ki,j “7/r , i 

soit en dérivant une nouvelle fois 

u M =u v* dans (Wi 


V i, j, k, l 


~'ij,kl 


■ £ — £ — £ 

kl,ij ik,jl jl,ik 


= 0 


ou encore 


Six équations 


^ £ 23,23 £ 33,22 £ 22,33 

£ 13,23 £ 32,31 _ £ 12,33 ~ £ 33,21 


+ permutation circulaire 
+ permutation circulaire 


(2.15) 


Réciproquement, si g vérifie (2.15), alors les formes différentielles 

dur.. = d: r, 

IJ I Kl, J JK,l I K 

sont exactes; elles permettent donc de construire le champ tu de tenseur antisymétrique. On vérifie ensuite 
que les formes différentielles 

du. = \tu., + dx, 

sont exactes, d'où la possibilité de construire un champ de déplacement u(X,t) défini dans iî ( . . 


3 Transport, dérivées particulaires 

3.1 Transport d'un volume 

Soit dü 0 un élément de volume de la configuration de référence, défini par trois vecteurs dX v dX 2 ,dX 3 . Par la 
transformation, ces trois vecteurs se transportent en trois vecteurs dx v dx 2 ,dx î qui définissent dans la 
configuration actuelle un volume dfl . 
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Figure 4 : Transport d'un élément de volume 

Le volume dfl est représenté par le produit mixte des vecteurs dx ] . dx 2 , dx 3 : 

dfl = (dx l A dx 2 j • dx 3 


donc 


dQ. — e.., dx , . dx„, dx 0 . 

ijk lj 2k oi 


Or, d'après (2.11) 

<Kl = e F. F, F dX, dX 0 dX„ 

ijk jp kq ir 1 p 2q or 

et, d'après (1.9) 

dÜ = e pqr det(F) dX lp dX 2q dX 3r = de^F)(dX 1 A dx\ ■ dX 3 

donc en définitive 


dn = Det(F) dÜ Q (216) 

3.2 Transport d'une surface orientée 

Soit dS un élément de surface de la configuration de référence de normale N . Par la transformation, cette 
surface se transporte en une surface ds de normale n dans la configuration actuelle. En considérant un 
vecteur V dans la configuration de référence qui se transporte en un vecteur v dans la configuration actuelle, 
on peut définir l'élément de volume ( dS N) ■ V qui se transporte en un élément de volume ( ds n)-v. 



D'après (2.16) 


ds n-v = det (F) dS N -V 


et comme avec (2.11) 

v = FV 


ds n ■ 

— 

FV 

= ds 

—T 

F n 






■7 = det F dS N-V 
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on obtient finalement 

dsn = det (F)F dS N 17) 

3.3 Dérivée particulaire d'une intégrale de volume 

Soit K(t) = ff] m k(x, t) dfl, une intégrale de volume sur le domaine Q(t) dans la configuration de référence. 

Pour en déterminer la dérivée temporelle, nous devons au préalable exprimer A'(f)sur la configuration de 
référence pour "passer" la dérivation sous l'intégrale. En effectuant le changement de variable (2.1), et en 
utilisant (2.16) 

dO = Det(F) dQ 0 = J dQ n 


on obtient 


m = fjf a j dn 0 


puis 


(ffil rrr 
dt - ^ 


dk T , dJ 
— J + k — 
dt dt 


dü„ 


A ce stade nous devons expliciter dJ / dt ■ En utilisant les notations indicielles, et en particulier les symboles de 
permutation, on a: 

J = det F = F F F 

g l ]k pqr ip jq kr 

soit 

dJ 1 9F 
= —e. e — —F. F 

dt 2 ,3k pqr dt n * 

or 


9F p d\{X,t) d 


dt dt dX dX 

P P 


du 


dt 


dX 




dv. dx, dv 

— — l - = —F 

dx. dX dx. p 

l p l 


donc 


mais 


dJ 1 dv. 

— = -£ e — FFF 

dt 2 * pqr dx, lp jq * 


£ FFF, = g det .F 

pqr lp jq lcr Ijk 


soit 


dJ_ 

dt 


1 dv. ~ 

— £..,£,., — L detF = 


ijk Ijk 


dx, 


6., — L det F 

il r\ 

dx, 


dv. 
— ^ J 
dx 

l 


dJ_ 

dt 


= J divv 


En reportant dans l'expression de dK / dt 


dK 

dt 



— J + k J divv 
dt 


dQ 


0 


(2.18) 
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puis en exprimant l'intégrale sur la configuration actuelle, on obtient finalement 


(f/t r r r 

dt ~~ s 


n (() 


dk , _ 

b k divv 

dt 


(2.19) 


En utilisant les égalités suivantes, 

dk dk 

— = b v-S/k 

dt dt 

div(kv) = v ■ V k + kdivv 

on peut écrire (2.19) sous la forme 


dt ~ s 


n(t) 


dk / , _\ 

b div kv 

dt K 1 


dü 


ou encore, en utilisant le théorème de la divergence 

= fffrn Jï dn + ff«w B H 


Application fondamentale: conservation de la masse 

La masse d'un système matériel qu'on suit dans son mouvement reste constante. 

M = JIJ m p(x,t)dn ^ = o 

et dt 

où p est la masse volumique. On a alors: 


— + pdivv — 0 ^- + div(pv) = 0 

dt ou dt V ’ 


( 2 . 20 ) 


3.4 Dérivée particulaire d'une intégrale de surface 

Soit K{t)= ff E(t) k(x,t)- n dS , une intégrale de volume sur le domaine E(t) dans la configuration de 
référence. Pour en déterminer la dérivée temporelle, nous devons au préalable exprimer K(t) sur la 
configuration de référence pour "passer" la dérivation sous l'intégrale. En effectuant le changement de variable 
(2.1), et en utilisant (2.17) 

dEn = det(F)F dT, g N 


on obtient 

K(t) = fJ Eo k(v(X,t),t)-J F dZ 0 N 


puis 


dK 

dt 


IL 


dk = ~ T - ^ 

d 

=~T 

' 

— ■JF N + k 


J F 

N 

dt 

dt 


, 


on doit donc calculer dF j dt 


F F — I 




dF 

dt 


F + F 



dF 

dt 



dt 
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donc 


dF ~ dv. dX dv . _ 

F = — 1 -e. <g>e = — L e.®e. =Vv 

dt dX dx. ' 1 dx. ’ 1 

k J 3 


dF 

dt 


F Vw 


= —V v F 


et 


dK _ „„ 
dt 


dk 1 - - =“ J _ - =- J _ 

— ■JF N + k- J divv F N -h- J X T v F N 
dt 




dK _ 


dk , . — » , , 

h divv k — « • V v 

dt 


J F Nd£„ 


puis en exprimant l'intégrale sur la configuration actuelle, on obtient finalement 


dK _ 
dt ~ ^ 


E(i) 


dk .. _ 7 __ 7 

h dîw k — vv k 

dt 


■ ndT, 


( 2 . 21 ) 


en utilisant la dérivée particulaire, (2.21) s'écrit 


dK r*r‘ 

Ht = JJm 


+ divv k + V k v —X7v k 
dt 


■ ndTi 


dF r- r- 

Ht = JJew 


+ ro< I fc A v I + v div k 
dt 


■ ndT, 


car 


rotik Av] = k divv — vdiv k + V k v — V v k 
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4 A retenir 

On appelle Variables de Lagrange le temps et la position initiale : X et t 
On appelle Variables d'Euler le temps et la position courante : J et t 


Dérivée particulaire 

Éâ = ÊA + (r.v\A 

dt dt \ 

Application linéaire tangente 

= dx . _ 

F — —e <g>e. 

dX ' 3 


Tenseur des déformations 

F (X,t)F(X,t)-l 


~ _ l 

~ 2 


Tenseur des déformations sous l'hypothèse des petites perturbations 


£ — — (v j w + Vu 


Transport d'un volume 

dFl = Det(F) dü 0 
Transport d'une surface 

dsn = det{F)F dS N 

Dérivée d'une intégrale de volume 


dK 

dt 


= JJI< 


sut) 


dk 

~dt 


div kv 


dQ 


Dérivée d'une intégrale de surface 


dK „„ 
dt ~ } 


m 


dk 7 

h divv k — V v k 

dt 


■ ndY, 
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EFFORTS DANS LES MILIEUX CONTINUS 


1 Définitions 

1.1 Forces 

Elles résument les effets mécaniques, autres que cinématiques, exercés sur le milieu continu considéré par le 
reste du domaine physique. Leur schématisation à chaque instant repose sur la définition d'un champ de 
vecteur $(x,t) et d'une mesure positive u> , définis sur la configuration actuelle f 2(f) . $(x,t) est une densité 
de force pour la mesure l o . 

* Si u) est une mesure de volume, alors $(J, t) est une force volumique (densité volumique de force) définie 
dans S! (t) de la configuration actuelle, par la fonction 

/: ÏÇfl(f)->/(ï,f)ë R 3 

* Si tu est une mesure de surface, alors $(f,£) est une force surfacique (densité surfacique de force) définie 
sur dü F (t) de la configuration actuelle, par la fonction 

F: x edn F (t)^ F(x,t)e R 3 

* ... etc ... 

Remarques: 

* Les forces sont définies sur la configuration actuelle. 

* A un instant donné et en un point donné te de d fi(f) , on ne peut imposer à la fois le déplacement et la force!. 
Mais l'un des deux doit être imposé. On note dQ F (t) la frontière où la force est imposée, et d£l F (t) la 

frontière où le déplacement est imposé. Dans le cas des appuis mobiles, les composantes non imposées 
cinématiquement le sont pour les forces 

* Le monde extérieur au milieu considéré doit, pour imposer le déplacement U(t) au bord dQ. (t) , exercer 

des forces que nous noterons R(x,t). Comme elles sont à priori inconnues, nous les appellerons réactions 
pour éviter de les confondre avec les autres forces qui, elles, sont données. 

1.2 Vecteur-contrainte et tenseur des contraintes 


1.2.1 Contrainte de Cauchy 

Soit un corps (C) en équilibre par application d'un système d'actions mécaniques 
extérieures. Imaginons qu'une surface E divise (C) en deux parties (1) et (2). La 
partie^,) est en équilibre sous les actions mécaniques extérieures qui lui sont 
appliquées et les actions mécaniques exercées par la partie (2). Nous admettrons que 
sur chaque élément de surface dT de E , (2) exerce sur (1) une force dF(x. t. rî) l/l} de 
densité superficielle T(x,t,n). 

dF(x, t, ti) 1/2 = f (x, t, n) dT, 



-37- 


Golay - Bonelli 




MMC 


f(x, t, n) est le vecteur contrainte au point x , relativement à la facette gE définie par son vecteur normal n . 
La densité surfacique de forces exercées en x dépend de x, t et aussi de l'orientation de la surface E au 
voisinage de x. Elle est linéairement dépendante de n . On introduit alors l'application a telle que: 

T(x,t,n) = o(x,t) n 

L'application a(x,t) s'appelle le tenseur des contraintes de Cauchy en i à l'instant t ; il caractérise, dans la 
configuration actuelle, les efforts intérieurs de cohésion exercés sur une partie du solide à travers l'élément de 
surface n dT, 


1.2.2 Autre écriture du tenseur des contraintes 
En utilisant (2.17), (3.1) devient: 

dF f x(X, t), t, n(N, t)] = n N(X) dS 


où II est le tenseur 

n( x, t): n e R 3 — *• n(x, t, n) = n(x, t)N e m 3 

défini par 

U(X, t) = (det F) o F (3 3) 

Cette application linéaire Ü(X, t ) , définie pour X e fl () , s'appelle le premier tenseur des contraintes de Piola- 
Kirchoff en X à l'instant t; la composante II est la \ ème composante du vecteur contrainte exercée sur la 
déformée d'une surface unité, normale à e., de la configuration de référence. On prendra garde au fait que le 
tenseur n n'est pas symétrique. 

Si maintenant on cherche le vecteur "force de cohésion" dans la configuration de référence 

dF 0 |x,f,ïvj = F (X, t ) dF p(X, t\ f n{N, t ) J = S N(X) dS 


où S est le tenseur défini par 


s = f n 


(3.4) 


Cette application linéaire S(X,t), définie pour X E fl 0 , s'appelle le second tenseur des contraintes de Piola- 
Kirchoff en X à l'instant t. Attention, sa composante S. n'est pas la i é ”' composante du vecteur contrainte 
exercée sur la déformée d'une surface unité, normale à e., de la configuration de référence, mais seulement la 
i eme composante de son transporté dans la configuration de référence. 

Selon le jeu d'écriture adopté, on a donc trois descriptions des contraintes: 


— 

— 

_1 — — = 7 1 

' — \ 

a — 

det F 

nf = 

det F 


F S F 


(3.5) 
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2 Equilibre 

2.1 Le Principe des Puissances Virtuelles (Germain 1972) 

Pour schématiser les efforts mis en jeu, il est commode d'imaginer des mouvements fictifs (ou virtuels) et 
d'analyser le travail ou la puissance qui en résulte. Par exemple, pour évaluer les forces de gravité agissant sur 
un objet, on peut imaginer de le soulever (mouvement virtuel de bas en haut). 

Un milieu matériel étant isolé, on peut distinguer les actions extérieures qui agissent sur le milieu, des actions 
intérieures qui représentent les liaisons existant entre toutes les parties du milieu. 

Axiome d'objectivité 

La puissance virtuelle des efforts intérieurs associée à tout mouvement rigidifiant est nulle. 

Axiome d'équilibre 

Pour tout milieu matériel repéré dans un référentiel absolu, à chaque instant et pour tout mouvement virtuel, la 
puissance virtuelle des quantités d'accélération n est égale à la somme des puissances virtuelles des efforts 

intérieurs n. et des efforts extérieurs n . 

2.2 Puissance virtuelle des efforts intérieurs 



Soit un milieu continu Q(t) d'intérieur O(t) et de frontière (t). Isolons maintenant un domaine £(f) de 
frontière <9£ (t) intérieur à Q (t), et soit n la normale en un point de d£(f) . A un instant f fixé, un mouvement 

virtuel défini par une vitesse virtuelle Sv est appliqué à £(f). Cette vitesse est supposée continue et 
continûment dérivable sur E(f) . 

Pour déterminer la puissance virtuelle des efforts intérieurs nous ferons les hypothèses suivantes: 

* IL admet une densité volumique p : 

11 -Jï.lyP dx 


* IL est en chaque point une forme linéaire des valeurs en ce point de dv et de ses dérivées premières: 


En décomposant le gradient des vitesses virtuelles en une partie symétrique SD et une partie antisymétrique 

SW, 


dëv 

dx 


SD + SW 


SW = - 

_ dSv 

SW. = - 

dSv. dëv . 

2 

dx dx 

<, 2 

dx dx. 

3 t ) 
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dëv 

dx 


dëv 

dx 



dëv. 

I 

dx 
l i 


dëv. 

i_ 

dx. 


la densité volumique des efforts intérieurs devient: 

p = A.ëV. + B. ëW.. — a ëD.. 

r l il IJ Jl IJ Jl 

Le premier axiome du principe des puissances virtuelles impose que pour tout mouvement de solide rigide la 
puissance des efforts intérieurs soit nulle. D'où: 

- Soit un mouvement de translation: ëv ^ 0 , èW = 0 et 6D = 0 


alors 


TT = fff'p dx = [I L A ■ ëv dx = 0 VE TT 
' JJJs ’ JJJe dans 

soit A ■ ëv = 0 V ëv , ou encore A = 0 

- Soit un mouvement de rotation: ëv = 0 , ëW ^ 0 et ëD = 0 
alors 

n,=/IU^= fff E B:6Wdx = 0 VE^q 
soit B : ëW = 0 V ëW , ou encore B = 0 . 


Donc en définitive: 

n , = ~ fff s ° ■■ 6D dx (3 6) 

On peut montrer que le tenseur a introduit ici correspond bien au tenseur des contraintes de Cauchy. 

2.3 Puissance virtuelle des efforts extérieurs 

Les efforts extérieurs comprennent 

- des efforts exercés à distance par des systèmes extérieurs à TT, supposés définis par une densité volumique 
de forces / , 

- des efforts de cohésion schématisés par une densité surfacique de force T sur dT, 

U a = fff, / ’ 6v dx + IL: T-ëvdx (3 y) 

2.4 Application du Principe des Puissances Virtuelles 

Si 7 est l'accélération et p la masse volumique de chacun des points de , alors 


n 


a 


dt 


fffsP v ■ ëv dx = fff z 


dpv 

~dt 


div(pv ® v ) 


• ëv dx 


n 


IL 


dv - dp ^ ■* ... 

p L v 1- p\7 v ■ v + div{pv)v 

dt dt 


■ ëv dx 


et en utilisant la conservation de la masse (2.20) et la définition de l'accélération (2.10) 

n „ = //Je P A ’ôv dx 


(3.8) 
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En application du Principe des Puissances Virtuelles on obtient: 

-//Je a :6D dx + JJL f ' 6v dx + fias T ■ 6v dx = fff E p 7 -6vdx (3 g) 

Pour exploiter le fait que (3.9) est vérifié pour tout mouvement virtuel, nous allons faire apparaître Sv dans 
chacun des termes. 

En appliquant le théorème de la divergence, le premier terme devient: 

-///e : 60 dx = -fff y a:^dx = -ff gs a-6v-ndx + [Jf y div x a ■ Sv dx 

Soit: 


//« 


T — a -n 


■6v dx + fffy 


f + diva - p 7 


•ôv dx V ôv 


Ou encore 


/ + div a = p7 dans E 
T = a ■ n sur <9E 


(3.10) 


2.5 Equilibre 


En considérant les développements du paragraphe précédent et en se ramenant au domaine fi(f), nous 
pouvons donc écrire les équations d'équilibre d'un solide soumis à un champ de forces extérieures / dansO(f) 

, à un champ de forces extérieures F e sur dQ F (t) et à un déplacement imposé IJ, sur dCl (t) ■ 

Dans la configuration actuelle: 


f(x, t) + div x a(x, t) = 0 V x E Q(t) 
a(x, t ) • n(x, t) = 


Fe(x,t) V x e di ï F (t) 

R{x,t) V x e dt \(t) 


(3.11) 

(3.12) 


Dans la configuration de référence: 


De même, si on note / 0 , R 0 et Fa les densités volumiques et surfaciques de forces mesurées dans la 
configuration de référence: 


f (X,t) + div x Il(X,i) =0 V x e fi g 


n (x,t)-N(x,t) = 


Fo(x,t) ViEi 1 (diî F (f), t) 

) VxedQ 0U 


(3.13) 


(3.14) 


Cas des petites perturbations 

Reprenons (3.10), en l'exprimant en fonction de X 

da 


f, [x(x, t), t ) + — n (x(x, t), t) = o v x{x, t ) e n(t) 

' ' ox 


, . da dx 

f (x(x, t ), t) + -^-(x,t) T -f- (X, t) = 0 v X e n 0 


dx. 
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Or x(X, t) = X + u(X, t) soit |f (X, t) = l + |f (X, t) 

oX oX 

On peut donc écrire l'équation d'équilibre sous la forme 


fi (x{x :, t ),t) + ^-(x,t) 


1 + -^{x,t) 

dX ’ 


= 0 


visa 


kj 


Sous l'hypothèse des petites perturbations, on peut alors écrire: 


d 'U / ~\r , . 

1 + (X, t) 

dX ’ 


âX 


soit 


, . da 


du. 

( x (X,t),t) + -^h(x,t) 

6 > 

~ ei « f > 

3 


Vieil, 


Enfin, en ne retenant que les termes d'ordre 0, et après avoir effectué un développement de f au voisinage de 
X, on obtient: 


da 


f,{ x ,t) + ^(x,t) = o Vieil. 


soit 


f(X,t) -f- div x a(X,t) =0 V x €E 


(3.15) 


Le raisonnement qui a permis de remplacer f\x(X,iyj par f(x,tj , permet aussi de remplacer F e (x(X,t)) par 
F e (X,t ) et R(x(X,t)) par R(X,t ) . Donc, comme condition sur la frontière on obtient: 


a{X,t)-N(X,t) 


Fe(X,t) 

R(X,t) 


v x e dn oF 
y xedn ou 


(3.16) 


2.6 Autre présentation: Principe fondamental de la dynamique 

(3.10) revient à écrire le Principe Fondamental de la dynamique. Dans un repère galiléen, pour tout système £ 

, le torseur dynamique (dérivée par rapport au temps du torseur cinématique) est égal à la somme des torseurs 
des actions intérieures. Soit: 


— f vdm = — îvpdX 
dtl J 


dt 

= f 

E 

= f 

E 

= f 


E 

dv p 


dt 


v pdivv 


dX 


dv - dp “• 

p F v 1- v pdivv 

dt dt 


P7 + v 


dp 

dt 


pdivv 


dX 

dX 


donc avec la conservation de la masse 
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— - f v dm — f pjdY, 

dt E e _ 

= //d£ + /a^<9£ 

E SE 

et le théorème de la divergence 

f p^dYi — f fdY, + J divadT, 

E E E 

on retrouve le bilan de la quantité de mouvement 

diva + / = P7 

L'équation de bilan sur les moments du principe fondamental de la dynamique s'écrit: 

///e 0MA ^ dm = //se OMAa-ndx + fff s OM A / dx 


(3.17) 


3 Quelques propriétés du tenseur des contraintes 

Dans tous les développements à venir, nous nous placerons dans le cas des petites perturbations pour un 
solide en équilibre. En conséquence, nous omettrons les variables set t. 

3.1 Symétrie du tenseur des contraintes 


On sait que 

///n 0MA (f J 7 ~7jdx = ff gn OM A an dx 

soit en notation indicielle 

IIL £ ijk x Mk ~fk)e,dx = ff gn e ijk x.a kl n l g, dx 

puis, par application du théorème de la divergence 


fff n 

fffn 


d 

~ A) ~ ^(VW e,dx= 0 

e # x jOr/ k ~f t ~ <V,) _ %%} e t dx=0 


et par application de l'équation du mouvement 

///nVV« dx=0 Vfi(f) 

c'est à dire 


e... a, . =0 Vz 

tjk kj 


ce qui implique 


£ 123°23 £ 132 Cr 32 ^ £ 213°13 ^ 231^31 ^ £ 312°12 £ 321 <J 2l ^ 

+ CT 23 - CT 32 = 0 - CT 13 + CT 31 = 0 + CT 12 “ CT 21 = 0 


donc en définitive 


a = a 

pq qp 


Le tenseur des contraintes est symétrique 
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3.2 Contrainte normale et contrainte tangentielle 



Considérons une facette de normale n . Tout naturellement, le vecteur contrainte 
T(n) peut être décomposé en une composante normale a n et une composante 
tangentielle r . 


et 


a = T(n ) ■ n — n ■ a ■ n 



( 3 . 18 ) 


( 3 . 19 ) 


On dira que cr est positive en traction et négative en compression. 

3.3 Directions principales, contraintes principales 

La matrice représentant le tenseur des contraintes est symétrique, elle est donc diagonalisable. Les valeurs 
propres sont réelles et appelées contraintes principales {cr v a n ,i t ui j . Les vecteurs propres, orthogonaux deux 

à deux, sont les directions principales ( nnnmnm ) ■ 0n a donc: 

<Ji=T{n I )-n I , a n =T{n n )-n n , a m = T(n m )-n m 


3.4 Invariants 


Le tenseur des contraintes possède trois invariants définis mathématiquement comme les coefficients de 


l'équation caractéristique det 


a — a 1 


. C'est à dire les quantité scalaires: 


= Tr(a) 


S " = 2 


Tr(af - Tria 2 ) 


= Det ( a ) 


Exprimés en fonction des contraintes principales, on obtient 


E r = °i + a n + a in 


S // - a i a n + a ii a in + a m a i 


S /// - a i a ii a iii 


( 3 . 20 ) 

( 3 . 21 ) 

( 3 . 22 ) 


3.5 Cercles de Mohr 

Connaissant le tenseur des contraintes a , on se propose de déterminer le domaine engendré par l'extrémité 
du vecteur contrainte quand n varie. Par commodité, nous nous plaçons dans une base orthonormée dirigée 
suivant les directions principales de a . Soit 
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n. 


G 

O 

b 


n a. 

i 

= 

i 

► 

i / 

n 2 

, o- = 

o 

o 

et T = 

U 2 a H 

TL 


b 

o 

g 


71, G 

3 


111 


3 III 


2 i 2 i 2 -i 

avec 7ij + n 2 + n 3 = 1 


on trouve aisément 


a 

n 


a, n; + a n n; + a m 


n 


2 

3 


et 


2 . 2 2 2 | 2 „2 | 2 2 
r + <? n = a I n x + n 2 + rc 3 

Dans l'hypothèse où les contraintes principales sont distinctes, on obtient alors après résolution du système: 


+ ( a n - g g )K - a m) 
0/ - 


r " + K - - a in) 

K - - a m) 


T ' 2 + K ~ g /)K ~ %) 
K/ - fT /)( f7 /// - %) 


Si on ordonne les contraintes principales de telle sorte que a i ^ a n ^ a „, , alors 


r 2 +(<r — u )(cr — a ) > 0 

t 2 +(<r — cr,) (a — a ) < 0 

r 2 + (<7 - (J )(a — a ) > 0 

'•n / ' ' n // ' — 


ou encore 


<7 — 


\ 2 

\ 

> 

a n “ a m 

2 


/ 


2 . 
T + 


< 


> 


(3.23) 


(3.24) 


(3.25) 



Dans le plan de Mohr, l'extrémité du vecteur contrainte, d'après (3.24), est donc 
intérieure au cercle centré sur Ou d'abscisse (c + a ) / 2 et de rayon 
(<7 — <7 ) / 2 . Par contre, d'après (3.23) (res. (3.25)), l'extrémité du vecteur 
contrainte est extérieure au cercle centré sur Oa d'abscisses (a„ + a riT ) / 2 
(resp.( (a l + a n ) / 2 ) et de rayon (a u - a m ) / 2 (resp. (a l + a n ) / 2 ). 
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Description des Cercles principaux: 

Nous allons étudier la description du grand Cercle de Mohr. Les facettes concernées 
sont parallèles à la direction associée à la contrainte principale a n . 

On constitue avec les directions I, III, Il un trièdre direct (0,e ,e ,e ) , la normale n 
de la facette évoluant dans le plan I III. 

Et on définit l'angle 0 = , et le vecteur t tel que (n,f ,11) soit direct. 

n = CosO e + SinO e 

et 

T = t 7 : CosO e : + a m SinO e m 

En utilisant les formules de changement de base de (0,e ,e ,e ) à (n,t ,11) , on a donc 



a i + 


+ 


Cos29 


Sin26 



Lorsque la facette tourne autour de la direction de la contrainte 
principale a d'un angle donné, l'extrémité du vecteur-contrainte 

tourne sur le cercle de Mohr d'un angle double dans le sens opposé 
(autour du centre du cercle). 
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4 Exemples de tenseur des contraintes 

4.1 Tenseur uniaxial 


a = a e 1 ® e 1 


a 0 0 
0 0 0 
0 0 0 


L'équilibre des forces sur la frontière du domaine nous donne: 



Sur E Q : n = —e l donc an — F o et Fa — —a e 3 


Sur E : : n = e 1 donc an — Fi et Fi = ae t 


Sur la frontière latérale les pressions sont nulles. 

On se trouve en présence d'un chargement uniaxial de traction/compression. 
Si a > 0 c'est un état de tension uniaxiale 
Si a < 0 c'est un état de compression uniaxiale 
La direction principale est e i 


4.2 Tenseur sphérique 


a = —pi 


—p 0 0 

0 —p 0 

0 0 —p 


Dans ce cas, toute direction est direction principale. La contrainte normale principale est -p. p est appelé la 
pression. Si p > 0 on a un état de compression, et si p < 0 on a un état de tension. 

Par exemple pour un fluide au repos: 

D'après l'équation d'équilibre 

diva + pg = 0 
—divpl + pg = 0 
—gradp + pg = 0 


soit 


dp 

dx 1 


dx 2 dx 3 


P9 


et P = P 0 ~ pgx 3 


Donc, pour un fluide au repos p + pgx 3 = Cste . 
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5 A retenir 


Vecteur contrainte et tenseur des contraintes de Cauchy 

T(x,t,n) = a(x,t ) n 

Le tenseur des contraintes est symétrique ! 

Equilibre 

f(x,t) + div x a(x,t ) = 0 Vie f2(t) 
u (x,t)- n(x,t) = 


F e (x,t) V x G dtl F (t) 

R(x,t) V x G dü [r (t) 


Contrainte normale 

a n = T(n ) ■ n — n ■ a ■ n 

Contrainte tangentielle 


J (— (-> — -.t 2 

a ■ n — n- a -n 
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ELASTICITE 


1 Approche expérimentale: essai de traction 

Pour déterminer l'évolution d'un système déformable, nous avons déjà déterminé les équations de la 
cinématique et de la sthénique. A ces équations, il est maintenant nécessaire d'adjoindre une relation 
supplémentaire reliant les efforts internes et les grandeurs cinématiques. Cette relation, appelée Loi de 
Comportement, dépend du matériau considéré. La construction d'une loi de comportement est basée sur des 
observations expérimentales. 

Dans ce chapitre nous exposerons le modèle de comportement des matériaux élastiques, sous l'hypothèse des 
petites perturbations. 




Pour effectuer un essai de traction simple sur un métal, on utilise une 
éprouvette cylindrique caractérisée par: 

- des extrémités surdimensionnées 

- des congés de raccordement (pour éviter les concentrations de contrainte) 

- une partie médiane cylindrique dans laquelle le champ de contrainte est 
supposé homogène, de traction simple parallèlement à l'axe de l'éprouvette. 

L'essai de traction consiste à enregistrer l'évolution de l'allongement relatif de la 
longueur initiale L g en fonction de la force de traction F, ou du rapport F / S 
, où S représente l'aire initiale de la section de l'éprouvette. 


Plasticité 



La figure ci-contre représente un tel enregistrement pour un acier 
inox. On remarque alors les propriétés suivantes: 

- Le diagramme est indépendant de la vitesse de chargement 

- La partie OA du diagramme est réversible. Si on charge jusqu'à un 
niveau inférieur à a 0 , alors la décharge décrit la même courbe OA. 

- La partie réversible est linéaire 

- Si on effectue un chargement au delà du seuil a , puis une 
décharge, l'éprouvette présente une déformation permanente. 

La partie réversible du diagramme de traction est, par définition, 
représentative du comportement élastique du matériau. cr 0 est la 

limite initiale d'élasticité du matériau. La linéarité du segment OA A 
caractérise le comportement élastique linéaire du matériau. 
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2 Loi de comportement élastique linéaire (en HPP) 


2.1 Forme générale 


A partir des observations expérimentales on peut écrire que les contraintes dépendent linéairement des 
déformations. En l'absence d'effets thermique et de contraintes initiales on a: 


a(x,t) = C(x) : e(x,t) 


(4.1) 


C est un tenseur du quatrième ordre, dont les composantes sont les coefficients d'élasticité du matériau. 

0 (•'•') -f- .. 

En utilisant les propriétés des tenseurs de contrainte et de déformation, on peut montrer que: 

ijkl ^ jikl ijlk ^ jilk 


Le tenseur C , dont la matrice représentative comporte 81 composantes, ne dépend donc plus que de 21 
paramètres indépendants. 

2.2 Matériau élastique homogène isotrope 


Toutes les directions sont équivalentes, de telle sorte que la loi de comportement est invariante dans toute 
rotation de la configuration de référence. Ce modèle s'applique à la plupart des matériaux: acier, béton, ... 


(4.2) 


Si la configuration est libre de contraintes, alors la loi de comportement s'écrit: 

<7 = A Tr(e) 1 + 2/jc 

ou encore en notation indicielle 

a. = Xe„6.. + 2 Lie.. 

ij kk i] ' ij 

Les coefficients matériel A et ^ , qui dépendent de la particule considérée, sont appelés les coefficients de 
Lamé. Leur expression en fonction du module d'Young E et du coefficient de Poisson v , est 


E 


ou 


A = 


E = 


2 (1 + v) 


et A = 


u E 


(1 + i/) (1 - 2v) 


/j,( 3A + 2 n) 


A + n 

avec, en inversant (4.2) 


et v = 


A 


e = - — Tr (°) 1 + 

h/ 


2(A + /ii) 


1 + v = 

a 

E 


(4.3) 


2.3 Matériau élastique homogène orthotrope 


Le matériau possède trois directions privilégiées deux à deux orthogonales. La loi de comportement est 
invariante par les symétries par rapport aux plans orthogonaux construits à partir de ces directions. Dans ces 
matériaux, on peut classer les tôles laminées, les composites tissés, le bois, certains bétons structurés, ... 

Dans ce cas on montre que la matrice de comportement est définie par 9 paramètres indépendants. Dans le 
repère principal d'orthotropie, la loi se met sous la forme: 
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1 

-"l2 

"13 

0 

0 



E , 

E, 

E, 


e n 


"21 

E, 

1 

X 

^23 

E 2 

0 

0 


^22 

^33 


“"si 

e 3 

-"32 

E 3 

1 

E 3 

0 

0 


2 £ 12 

2^3 


0 

0 

0 

1 

G 12 

0 


2^3 


0 

0 

0 

0 

1 

G 23 



0 

0 

0 

0 

0 

les conditions de symétrie 




"12 

" 21 

"13 

"31 

V V 

32 23 



E , 

E 2 

E , 

e 3 

E 3 e 2 




G 


(4.4) 


2.4 Matériau élastique homogène isotrope transverse 

Un matériau homogène isotrope transverse est tel que la matrice de comportement est invariante par toute 
rotation autour d'un axe privilégié. En utilisant cette invariance, on montre que seuls 5 paramètres 
indépendants caractérisent le comportement. Si l'axe est porté par la direction 3, on a alors: 


2e, . 


1 

-"12 

-"13 

E x 

E r 

E i 

^21 

1 

-"13 

E 1 

E r 

E i 

-"31 

-"31 

1 

E., 

E., 

E, 

3 

3 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 
0 

— 0 


G, 


12 

0 


G, 


G,. 


(4.5) 


2.5 Caractéristiques de quelques matériaux 

Matériaux isotropes usuels: 


Matériau 

E en Gpa 

V 

p en kg/l 

acier 

210 

0.285 

7.8 

fonte grise 

90 à 120 

0.29 

7.1 

aluminium 

71 

0.34 

2.6 

béton 

10 

0.15 

2.4 

fibre de verre E 

73 

0.15 

2.54 

Graphite HM 

350 

0.4 

1.92 

résine époxy 

3.8 

0.31 

1.15 
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Matériaux composites: 



Unidirectionnel 

Tissu 

Unidirectionnel 

Unidirectionnel 


Verre/Epoxy 

Verre/Epoxy 

CarboneHT/Epoxy 

Kevlar/Epoxy 


50% 

50% 

50% 

50% 

p en g/cm 3 

1,87 

1,87 

1,49 

1,32 

E l en Mpa 

38000 

21000 

116000 

65000 

E , en Mpa 

11500 

21000 

7500 

4900 


0,28 

0,26 

0,32 

0,34 


2.6 Critères de limite d'élasticité 

Les critères de résistance que nous allons définir représentent des valeurs limites pour les contraintes 
maximales, et permettent de ce fait de garder un caractère élastique aux déformations. 


2.6.1 Critère de Tresca 

Il consiste à considérer de manière indépendante les trois contraintes de cisaillement maximal du tricercle de 
Mohr. Soit en fonction des contraintes principales 


Sup 


ny 


' ni ’ ru 


< 2 <7„ 


(4.6) 


2.6.2 Critère de Von-Mises 


77 K a / - a nf + K - G iJ + ( a n - a m) 2 


< 


ou encore 



2.6.3 Critère de Hill 

le critère de Hill s'applique dans le cas de matériaux élastique orthotropes 

F ( a n - + H ( a n - a 33 f + G K - a 3 J + 2 La l 3 + 2Ma n + 2N = 1 

où F,H,G,L,M,N sont des constantes fonctions des contraintes à ruptures. 


(4.7) 


(4.8) 
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3 Le problème d'élasticité 

3.1 Ecriture générale 

Cinématique : + Equations de compatibilité 

£ = iJvu + V r w) 

u = Uo(X) sur 

Equilibre : 

div a + / = 0 dans Q 


F 


an = 


R 


sur dfl p 
sur dfl i: 


Loi de comportement : 

a = A Tr(e)I + 2/ie 

3.2 Formulation en déplacement 

diva + f = 0 

div (A Tr{e)I + 2 fie) + f = 0 

A V(Tr(e)) + 2 fjdiv(e) + f = 0 

A S/ [div u) + ndiv{y u) + fidiv(V J u) + f = 0 


soit l'équation de Navier 

(A + jt i)S7(div u) + fidivÇV u) + f = 0 


Remarque : Si on prend la divergence de l'équation de Navier 

(A + 2/j,)A(div u) + div(f) = 0 


(4.9) 


Donc, si le champ de forces volumiques est tel que divf = 0 alors div u est une fonction harmonique. 


3.3 Formulation en contrainte 

En partant de l'écriture des équations de compatibilité, on peut démontrer les équations de Michell 

divÇ V a) + — — V( V Tr(a) ) + — — divf / + V/ + V r / = 0 
1 + v 1 — v 

Soit, si le champ de force est uniforme, on obtient les équations de Beltrami. 

(1 + v)div(X7 a) + V(V Tr(a )) = 0 

3.4 Théorème de superposition 
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Si ( U,f,F ) et (V, g, G) sont deux jeux de données engendrant respectivement des solutions u et v , alors 
au + /3v est solution du problème de données (aU + /3V, af + f5g, aF + (3G) (Le problème est 
évidemment linéaire). 

3.5 Elasticité plane 

3.5.1 Contraintes planes 

Dans le cas où le chargement est dans le plan 12, la structure mince dans la direction 3, on peut faire 
l'hypothèse que le problème est plan et libre de contraintes dans la direction 3. 

Dans ce cas 


On remarquera que la déformation suivant 3 est non nulle. 

3.5.2 Déformations planes 

Dans le cas où le chargement est dans le plan 12, la structure très élancée dans la direction 3, sans possibilités 
de déplacement suivant 3, on peut faire l'hypothèse que le problème est plan sous l'hypothèse des 
déformations planes. 

Dans ce cas 


= ct h(^z 2 ) cr 12 (x v x 2 ) 0 
<j = o 21 (x v x 2 ) cr 22 (x v x 2 ) 0 


0 0 0 


et d'après la loi de comportement 


= e n( x v x s) e n( a 
e= ejx^xj e 22 (a 



0 


= s n( x v x 2) e i2 (*!>**) 0 

£= £ 21 (x v x 2 ) £ 22 (x v X 2 ) 0 


0 0 0 


et d'après la loi de comportement 


= a il( X V X 2) 

( 7 = crjx^xj <t 22 (c 



0 


On remarquera que la contrainte suivant 3 est non nulle. 
D'après (4.3) 



donc 


33 



Nous allons prouver que les contraintes peuvent être déterminées par une seule fonction scalaire. 
En appliquant l'équation d'équilibre (3.11) on a : 
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Vi + Œ V. 


= o 
= 0 


donc 


3 0(*i>* 2 ) /a n =cj) 2 et cr 12 = -^ 
3 ^{x v x 2 ) / 0' 21 ='0 2 et cr 2 = 


comme le tenseur des contraintes est symétrique, ona^+^-O, donc 


^X(x v x 2 ) /<t> = Xo et -0 = 

en définitive on a prouvé 


-X,i 


3x0^ S ) / 


a n = X 22 
C7 22=X 11 

CT 12=-X 12 


X est appelée la fonction d'Airy. 

Le tenseur des contraintes devant vérifier l'équation de Beltrami (4.11), on a 

( 1 + ^K«+Vv = 0 

d'où 

AAx = 0 

X est donc une fonction biharmonique. 

3.6 Thermoélasticité 


3.6.1 Thermodynamique : équations de bilan 

Jusqu'à présent nous avons utilisé les équations de bilan suivantes: 

Conservation de la masse 

dp -* 

b p divv = 0 

dt 

Conservation de la quantité de mouvement 

div a + f = pj dans SI 

Conservation du moment cinétique (3.17) 

Nous introduisons maintenant l'équation de bilan de conservation d'énergie, ou encore le premier principe de 
la thermodynamique: 

f( £+A ') = p - + « 

où 

E représente l'énergie interne E — f pe dfl (e densité d' énergie interne) 

n 

K représente l'énergie cinétique K = f^pv-vdiï (v la vitesse) 

O 
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P représente la puissance des efforts extérieurs P t = f f -v dfl + f F ■ v dfl 


an 


Q représente le taux de chaleur reçu Q = J r dfl — f q - n dfl (g vecteur de chaleur et rsource de chaleur) 

n an 

Par application du premier principe, en utilisant (5.7) on a: 

(j p 

f p — dfl + f pv ■ 7 dPl = J f ■ v (Kî+ J F -v dQ + JrdQ — J q-ndQ 

n dt n n an n an 

en utilisant la conservation de la quantité de mouvement (3.9) 

fj p 

J p - — dfl — f a : e dfl + f r dfl — f q ■ n dfl 

n dt n n an 

Soit, par application du théorème de la divergence, la forme locale du premier principe 

pè = a : £ + r — divq 

Nous présentons également, sans plus de discussion le second principe de la thermodynamique: 

,L iü .,î£l dü 

dt n T an T 

où T est la température et S l'entropie. Ce second principe s'écrit sous sa forme locale 

q r 

ps + div >0 

T T ~ 

où s représente l'entropie massique 


(4.12) 


(4.13) 


3.6.2 Equation de la chaleur 

On peut exprimer l'énergie interne massique e en fonction de l'entropie massique s, de la température T et 
l'énergie libre tp . 


e = ip + Ts 


(4.14) 


En thermoélasticité, sous l'hypothèse des petites perturbations, pour un écart de température par rapport à la 
température au repos T _T o 

petit, on a: 


-0 = T) 

Grâce au second principe on peut montrer que 

~ dip dtp 

(j = p — — et s — 

d~e dT 

donc, le premier principe peut s'écrire 

pè = pip + pTs + pTè 

et comme 


; d^P = 

ip = : e ■ 

de 


— = — : é — sT 
dT p 


on a 


a : e — psT + psT + psT — a : ê + r — divq 


or 
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dip' 

dT 


d 2 i/> ,7 d 2 ip ÿ _ 1 do 7 | Os ÿ 

9 = 9T - 5T 


c'est à dire que le premier principe s'écrit 


-T — - : è + pT — - T — r — divq 
dT dT 

En introduisant la chaleur spécifique C = T 

OT 


—T —— : è + pCT — r — divq 
dT 


puis la loi de Fourier q = — kS7 T , où k représente la conductivité thermique, 

O “ “ 

-T— :ê + pCT = r + divkV T 
dT 

En général la contribution mécanique est négligeable par rapport aux autres contributions, si bien que 
l'équation de bilan de l'énergie conduit à l'équation de la chaleur : 


pCT = pC 


dT 

~dt 


+ v-VT 


— r + divkV T 


(4.15) 


Dans le cas où le problème à traiter est stationnaire, sans source de chaleur, avec une conductivité constante, 
on retrouve l'équation habituelle : 

AT — 0 


3.6.3 Loi de comportement thermo-élastique 

Dans le cadre de la thermoélasticité , l'énergie libre spécifique s'écrit comme un développement limité au 
second ordre en déformation et température, ou plutôt en déformation et écart de température t = T — T a 
(supposés " petits ") : 


M Ê > T) = ^e : C :s- psr-^pbr 


(5 : £T 


Par définition 


a = p^t(s,T) = C : e- (3 t = C : 
de 



où a représente le tenseur des dilatations thermiques 

Dans le cas isotrope la loi de comportement thermo-élastique s'écrit : 

a = \Tr(e)l + 2 pe — (3A + 2 p)ocr 
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4 A retenir 


Loi de comportement élastique linéaire isotrope 

<7 = A Tr(e) 1 + 2 p e 


v — - 1 ~\~ V — 

e = Trio ) 1 + a 

E v E 


Critère de Tresca 

Sup 


< 2 <7„ 


Le problème d'élasticité 


1 


e = ^-|V u + V 1 u 
u = Uo(X) sur dü v 
div a + / = 0 dons il 
F sur dfl„ 


a n = 


R sur <9L1, 


cr = A Tr(e)I + 2pe 


Equation de Navier 

(A + p)V(divu) + pdiv(y u) + f — 0 


En élasticité plane sous l'hypothèse des deformations planes : 


CT 33 = V K + CT 22 | et £ 33 = 0 


Conservation de l'énergie 


d 


-( E + K) = P„ + Q 

Forme locale de la conservation de l'énergie 

pè — a : e + r — divq 

Equation de la chaleur 


P C 


dT - 

+ v-VT 

dt 


— r + divkV T 


Loi de comportement thermoélastique isotrope 

a = XTr(e)l + 2ps — (3A + 2 p)ar 
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INTRODUCTION A LA MECANIQUE DES FLUIDES 


1 Loi de comportement 

En mécanique des fluides, nous travaillerons toujours en variables d'Euler 

Comme pour les matériaux solides (qui sont des fluides qui s'ignorent ..) les lois de comportement fluide sont 
élaborées à partir de l'expérience. 



1.1 Fluide Newtonien 

Pour un fluide Newtonien, les contraintes sont une fonction affine des vitesses de déformation. Soit, 

a = —pi + XTr(D)I + 2 pD 


ou 


D = -\gradv + grad 


soit, en notation indicielle 


1 

dv. 

dv. 
_| L 

2 

dx. 

J 

dx. 


H est la viscosité dynamique (dimension Poiseuille = -M-) 

LT ' 

A est le second coefficient de viscosité 


On introduit également la viscosité cinématique v = — (dimension Stokes= — ) 

P T 


._ L 


(5.1) 


(5.2) 
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1.2 Fluide incompressible 

Si le fluide est incompressible, alors on a vu que divv = 0 ou TrD = 0 


Donc, (5.1) devient 


a = —pi + 2fiD 

1.3 Fluide non-visqueux 

Si le fluide est parfait, alors on a ne tient pas compte de la viscosité, donc (5.1) devient 


Le tenseur des contraintes est alors sphérique. 

En particulier, l'action d'un fluide non visqueux sur une paroi est normale à la paroi (d'après l'équation 
d'équilibre). 

1.4 Fluide au repos 

Si le fluide est au repos, alors v = 0 , donc (5.1) devient 


a — —pi 


(5.5) 


2 Conservation de la masse 

La masse d'un système matériel qu'on suit dans son mouvement reste constante. 

M = film P& l ) dÇl et = 0 

où p est la masse volumique. On a alors (2.20) 

— + pdivv = 0 ou — + div(pv\ = 0 

dt dt \ I (5-6) 

Si on considère une grandeur différentiable '5 quelconque, on a alors pour un fluide incompressible 

j t fff m 1 Pdrn = JJf m ^dm {5 7) 

Démonstration: 


J t fff m Tp dm = Jffn(t) ^ P dÇl 


= fffm 

= fffn » 

- JIU 

= fffm 


dipp 


dt 

dit 

o — 

dt 


ippdiv v 
dp 






dt 

dp 

dt 


dip 

dt 


dit 


fffsi(t)~ dm 

jjjm dt 


dfl 

xppdiv v 
pdivv 


dfl 

dfl 


Soit E un domaine géométrique fixe traversé par le fluide. 
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fias P v -nddZ = fff E div(pv) dZ 

= fff s -—dZ 

JJJz dt 




(i d'après le théorème de la divergence ) 
(d'après la conservation de la masse ) 

(car Z est fixé ) 


Si le fluide est incompressible, alors la masse volumique est constante et 

dp_ = dp_ = 0 
dt dt 

Si on note g le débit massique à travers une surface Se t g le débit volumique, alors 

Ç m = ff s pv-nddZ = pff s vnddZ = pq v 

Donc, en définitive: 

Pour un domaine Z fixe traversé par un fluide incompressible ff. )r v-n ddZ — 0 : le débit volumique à 
travers la frontière dZ est nul. 


3 Equation du mouvement 

D'après l'équation du mouvement(3.10), 

7 = dv 

t + diva — p — 
dt 

D'où, pour un fluide newtonien 

dv 


dt 


f + div 


-pi + A Tr(D)I + 2 pD 


= f — div 


pi 


Xdiv 


Tr(D)I 


2pdiv 


D 


/ - Vp + AV 


fdiîwj + pdiv | 


Vv + V 1 v 


Soit l'équation de Navier-Stokes compressible 


dv 




p — = f — Vp + (A + /i)V [divv] + pAv 
dt \ I 


* Pour un fluide incompressible , divv = 0 , donc (5.8) devient: 

= f — Vp + pAv 


dv 

p Tt =p 


dv - ^ 

h V ■ VU 

dt 


* Pour un fluide non visqueux, (5.8) devient: 


dv 

p Tt =p 


dv - ^ 
dt. 


= /-Vp 


(5.8) 


(5.9) 


(5.10) 


-61- 


Golay - Bonelli 



MMC 


4 A retenir 

Loi de comportement pour un fluide newtonien 

a = —pi + XTr(D)I + 2 pD 

Conservation de la masse pour un fluide incompressible 

divv = 0 

Grâce à la conservation de la masse pour un fluide incompressible 



Equation de Navier Stokes compressible 


p — — f-Vp + (X + p)V + pAv 


dt 
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Annexes 


ANNEXES: RAPPELS DE MECANIQUES DES SOLIDES RIGIDES 


1 Cinématiques du solide 

Avertissement: L'objectif de ce chapitre, est de familiariser les étudiants avec les notations tensorielles. Afin 
d'en simplifier le contenu, nous ne considérerons que des bases orthonormées. 

1.1 Description du mouvement 

Soit S un ensemble de particules tel que la distances entre deux particules quelconques reste pratiquement 
constante au cours du mouvement. On étudie l'ensemble S en le considérant indéformable: solide rigide. 


1.1.1 Système de référence 

Dans un espace euclidien £ à trois dimensions, soit e , e 2 , e 3 une base orthonormée. On définit un référentiel 

d'observation par cet espace euclidien et le temps: . On définit la dérivée d'un vecteur U par rapport 

au temps dans ce référentiel par: 


dU dU. _ 

dt dt ' 


1.1.2 Mouvement d'un solide 

Soit S un solide rigide en mouvement par rapport à K . Soit Ç s (o,e v e 2 ,e^j un espace euclidien lié à S. 

Considérons un vecteur lié à S, dont les composantes sont représentées par X dans Ç et Is dans £ s . 
On note Q l'opérateur linéaire définissant le passage de £ à . 

X s = QX et X = Q t X s 

Comme Xs est indépendant du temps puisque lié à S, on a: 


dX 

dt 


dQ T 

dt 


Xs = 


dQ T 

dt 


QX = LX 


or 


Q T Q = I 


C'est à dire 


^-Q + Q T ^-= 0 
dt _ dt 

L + L t = 0 

L est un opérateur linéaire antisymétrique, on peut donc définir un vecteur O (d'après (1.10)) tel que: 


dX 

dt 


= LX = fl A X 


( 6 . 1 ) 


avec 
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fl = —e...L..e. 

2 y* 3 i k 


1.1.3 Torseur cinématique 

Soient A et P deux particules du solide S. 

OP = ÔA + ~ÂP 

donc par dérivation 


dOP 


dt 


dOA 


dt 


+ 


dAP 


dt 


dOA 


dt 


+ n a ap 


soit 


v(p) = v{a)+üaap 

On définit le torseur cinématique par le vecteur vitesse de A par rapport à K , V (j4) et le vecteur de 
instantanée fl : v| = 


, 

V[A) 


n 


1.1.4 Accélération 

V(P) = v{A)+flAJp 


Soit 


dv[p ) 

dV(A) 

dfl 

dt 

dt 

dt 





-, s 4Q 


->( F ) = 

7 U)+^- 

A AP 

v ' dt 



A A P -f Q A 


dAP 


dt 


A AP + n A A AP J = 7(^4 

1.2 Composition des mouvements 


dfl 

dt 


A AP + \ fl ■ AP\fl - fl 


L'espace temps est commun à tous les référentiels d'observation, on considère deux référentiels üR“ 

5R 6 (Ç ,i). 


1.2.1 Dérivation composée 

Soit U un vecteur dans la base (0,e 1 °,e“,e“j , U = U.e“ par dérivation: 


dU 


dt 





Soit U un vecteur dans la base £ 


s~\ -*b ~*b ~*b 

O, e 1 , e 2 , e 3 


U = U £ 


par dérivation: 


dU dU _ b TT de* 

= L e. + U , — - 

dt dt ' ' dt 

sr 


( 6 . 2 ) 

rotation 


AP 

(6.3) 


C,t) et 
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Car Ç b est en mouvement par rapport à , et d'après (6.1) 


de 77 

— — = Sleve 1 A e b 
dt ' 

d'où, avec fi^/e 4 représentant le vecteur rotation de Ç 4 par rapport à £“ 

+ lV/f‘ A U 


dU 

dt 


dU 

dt 


1.2.2 Composition des vitesses 

Soit P £ S : 


O a P = O a O b + O b P 


dO“P 

dO a O b 

dO b P 

dt 

dt 

3î a 

dt 

3Î“ 


Et donc d'après (6.4) 


soit 


sr 


dO“P 

dO“O b 


dO b P 

dt 

dt 

5ft D 

dt 


+ n?/e a o b p 


Va(P) = Fi (P) + Ve (P) 

Vitesse /3?“ Vitesse /5R b Vitesse d ’ entrainement 


(6.4) 


(6.5) 


1.2.3 Composition des accélérations 
On dérive (6.5) par rapport à : 


d 2 O a P 


dt 2 

» 

la{P) = 


d O a O 
dt 2 


d 2 O b P 
dt 2 


+ A 


dO b P 

dt 


dQ? /f l 

dt 


a o b p + a 


dO b P 

dt 




ia(o b ) + j b (P) + ^/e AV b (P) + 


/ f 

/ç‘ 

+ fV /ç l A /ç‘ 

A O b P + -(-fV/e* A 

d0 4 P 

+ fV/î* A 0 4 P 

dt 

dt 

s 1 


K* 


iSP) = i b {P) + iAo b ) + 


dQç>/ç b 


dt 


A 0‘P + fievï* A /{* AO‘P +2f2 ?V f‘ A ^ (P) 


Accélération de Coriolis 


Accélération d ' entrainement 


7u(P) = 7i(P) + 7e(P) + 7c(P) 


( 6 . 6 ) 
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2 Cinétique 


La cinématique ne s'intéresse au mouvement des corps que du point de vue de l'espace et du temps: durée, 
vitesse, distance, etc ...; tandis qu'en cinétique on introduit, en plus, le concept de masse, c'est à dire qu'on 
tient compte aussi de la masse 

2.1 Définitions 

On définit la masse d'un solide S par : 

m = fffs dm ( P ) = fffs P( P ’ L ) dv (6. 7) 

Où p représente la densité volumique de masse. 

On définit G le centre de masse (ou d'inertie) du solide S par : 

v O e £ fff s OPdm(P) = mOG (g g) 

2.2 Eléments de cinétique 


2.2.1 Torseur cinétique 

On définit le Torseur Cinétique ou Torseur des quantités de mouvement par : 

R = Jff s V(P)dm(P ) Résultante cinétique de S /St 

kA = fffgAP AV(P)dm(P) Moment cinétique en A /St 

On peut remarquer que si le repère est fixe, alors : 

R = ~ fff s OP(P)dm(P) = mV(G) 

2.2.2 Torseur dynamique 

On définit le Torseur dynamique par : 

d = JJJ j(P)dm(P) Résultante dynamique de S /St 
Sa = fffgAP A 7 (P)dm(P) Moment dynamique en A /St 




( 6 . 9 ) 


( 6 . 10 ) 


2.2.3 Relation entre torseur cinématique et torseur dynamique 
En dérivant par rapport au temps dans K on obtient : 


2_dR 

dt 


( 6 . 11 ) 


^ = fffs A V(P)dm(P) + fffg AP A j(P)dm(P) 

= fffs A HP)dm{P) + fffg ^ A V(P)dm(P) + fff s AP A j(P)dm(P) 

= -V(A)Afff s V(P)dm(P) + ÔA 


Sa = — + V(A) A mV(G) 
dt 


( 6 . 12 ) 
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Sa = M s AP A j(P)dm(P) = fff s AG A 7 (P)dm(P) + fff s GP A 7 {P)dm(P) 
Sa = AG A d + 6 g 


Sa = — + AGAm 7 (G) 
dt 


( 6 . 13 ) 


2.2.4 Energie cinétique 

On définit l'énergie cinétique du solide S par : 

T(S) = ±fff s V(P)dm(P) 


( 6 . 14 ) 


2.2.5 Théorème de Koenig 

Soit ^(0,e v e 2 ,e 3 ^ un espace euclidien et un espace euclidien barycentrique lié au solide S. 


k A = ko + AG AmV{G) 
Sa = 6g + AG A m 7 (G) 


( 6 . 15 ) 

( 6 . 16 ) 


T(S) = fffs — -— dm 
’ JJJs dt dt 


1 ,,, dOG dOG , 1 ,,, dGP dGP , ,,, dGP dOG , 

T(S) — — fff a dm H — fff dm P fff. dm 

2 JJJs dt dt 2 JJJs dt dt JJJs dt dt 


1 


T (S) = Tjj (S) + -mV 2 (G) 


2.3 Cinétique du solide rigide 


( 6 . 17 ) 


2.3.1 Opérateur d'inertie 
On définit l'opérateur d'inertie par J A tel que : 

3 A \ u e £ — > J A (u) — fff AP a(u A APjdm 


Si AP — xe alors 


AP A u A AP = £...x.£. u x e. = S .6 x u x e. — S .6 x u x e. = x.u.x.e. — x.u.x.e 

l / ijk j kpq p q i pt qj j p q t qi pj j p q t J i J i J J i i 

AP A (u A AP^i = (are, ® e,) • («.e.) - (x k x.e k ® e.) • (u.e.) = (x 2 .e k ® e, - x k x.e k ® e.) • (^e ) 

Et l'opérateur d'inertie est représenté par la matrice : 


h = 


h -Aa 

12 z _ ^23 

_ ^13 _ ^23 3 


OÙ 
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I l=fff( X 2+ X 3) dm 1 12 
h = IIl( X l+ X l) dm J 13 
d 3 = fff( X l +^t] dm d 23 


///* l X 2 dm 

fff x i x 3 dm 

fffx 2 x 3 d m 


2.3.2 Influence des symétries matérielles 

• Si le solide S possède un plan de symétrie l^A,e v eA , alors 

fff X l X 3 dm = fff X l X 3 dm + fff X l X 3 dm = fff X l X 3 dm ~ fff X l X 3 dm = 0 


Soit 


h = 


h - J 12 « 

h « 

0 0 7 


Si le solide S possède un axe de symétrie ( A, e 3 ) , alors 


Et comme 


///••• pdx t dx 2 dx :i = fff- prdrd,0dx 3 

fff x 1 x 3 pdx 1 dx 2 dx 3 = fff r cos 9x 3 prdrd9dx 3 = 0 

si a 

fff x 1 x 2 pdx 1 dx 2 dx 3 = fff r 2 sin 9 cos 0x 3 prdrd0dx 3 = 0 


on a finalement 


h = 


f 0 0 

0 I 2 0 

0 0 7 , 


• Moment d'inertie par rapport à une droite A (de vecteur unitaire S ) passant par A 
Soit H le projeté orthogonal d'un point P du solide S, on a alors : 


7 a = fff PH dm = fff 


AP 


AH 


dm = fff 


AP 




dm 


h=fff 


6.6 \ AP. AP \- AP. 6 


dm = fff 6. ~s\apAP^-Ap{aP£ 


dm 




Et comme 


= d -fff d^AP.Apj-AP^AP.^ dm 
a A (b A c j = (a ■ c j b — (d ■ b j c , 

= 6. fff ÂP A p A Âpj dm = 6 ■ J A {6) = 6 ■ 7 , • S 


d .= 
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Soit 


soit 


Théorème de Huyggens généralisé 

J A {u)= fff AP A (m A APj dm 
= fff AG A A Apj dm + fff GP A A Apj dm 
= A~G A A fffÂPdmJ + fff GP A (u A Ig) dm + fff GP A (u A GP^dm 
= ~ÂG A |u A mlc) + (fffGPdmj A (u A ÂgJ + fff GP A |u A Gpj dm 


J A (Ü) = AG A\u A mA GJ + J G (Ü) 

Théorème de Huyggens appliqué au moment d'inertie par rapport à une droite A (de vecteur unitaire 
6 ) passant par A (tel que A G _L <5 ) et A' (de vecteur unitaire 6 (passant par G. 


L =6- J.(s) = srs = s- 


AG A\ë AmAG\ + J J6) 


= 6 - 


mAG 6 


-miô-AG^AG + lJ 


^ = SI A S = S- 


mAG 


6 = mA G + I„ 


2.3.3 Moment cinétique du solide 

L = fff g AP A V(P)dm(P) 

Soit Q un point quelconque du solide 

L = fff s AP A [v(Q) + ÜAQP^ dm(P) 

Ia = mAG A V(Q) + fff g {Jq + QP^j A A Qpj dm{P) 

k A = mÂGAV{Q) + ÂQA^A fff s QPdm(P)j + fff s QPA^2A Qpj dm(P) 

D'où 

L = mAG A V{Q) + mAQ A (fi A Qg) + J G (Ü) 

Dans le cas particulier où Q=G, on obtient : 

L = mAG AV(G) + J G (fi) 

C'est-à-dire 

ko =J a (i ï) 

2.3.4 Energie cinétique du solide 

T(S) = ±fff g V(P)dm(P) 

Soit Q un point quelconque du solide 
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Soit 


T(S) = - fff s V(Q) + fiA QP) dm(P) 


T (S) = -ffI s v (Q)j dm(P) + ffflv(Q).ÜAQPdm(P) + -fff s ÜAQP dm(P) 


T (S) = — V (Q) + mV(Q)-tSlAQGj + -fff s Sl 


QP a\ÜAQP 


dm(P) 


T (S) = (Q) + mV{Q) ■ [fi A QGj + ■ J q (Q) 


Dans le cas particulier où Q=G, on obtient : 


T(S) = jV(G) + ^.J G (à) 


3 Equations fondamentales de la mécanique des solides 


3.1 Torseur associé aux efforts externes 

Soit /(P) une densité volumique de force exercée sur le solide S . 

Soit F (P) une densité surfacique de force exercée sur la frontière du solide dS . 
Soit #(P) une densité linéique de force exercée sur une courbe T . 


Soit Fi une force ponctuelle exercée en un point P de S . 


Le torseur des efforts extérieurs est défini par : 


US) 


r = III m + il m + 1 m + z f, 

s + as r * 

Ca = fffAPAf(P) + ffAPAF(P) + fAPA5(P) + ZAPAFi 

s as r * 


3.2 Loi fondamentale de la dynamique 

Il existe au moins un référentiel Galiléen associé à une chronologie, tel que : 

VjS, Vf Torseur dynamique — Torseur des forces extérieures 


Ou encore 

VS, Vf Vl(d,^) = 5 e (P,Cx) 

En conséquence, on peut énoncer : 

Théorème de la résultante dynamique : dans un référentiel galiléen 

R = m"){G) 

Théorème du moment dynamique : dans un référentiel galiléen, soit A un point fixe 

dk A -pi 

oa — — La 

dt 
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